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We study an inflationary scenario with a two-form fleld to which an inflaton couples non-trivially. 
First, we show that anisotropic inflation can be realized as an attractor solution and that the two- 
form hair remains during inflation. A statistical anisotropy can be developed because of a cumulative 
anisotropic interaction induced by the background two-form field. The power spectrum of curvature 
perturbations has a prolate-type anisotropy, in contrast to the vector models having an oblate- 
type anisotropy. We also evaluate the bispectrum and trispectrum of curvature perturbations by 
employing the in-in formalism based on the interacting Hamiltonians. We find that the non-linear 
estimators /nl and tnl are correlated with the amplitude of the statistical anisotropy in the power 
spectrum. Unlike the vector models, both /nl and tnl vanish in the squeezed limit. However, the 
estimator /nl can reach the order of 10 in the equilateral and enfolded limits. These results are 
consistent with the latest bounds on /nl constrained by Planck. 



I. INTRODUCTION 



■ The inflationary paradigm [l| can successfully account for the observed temperature fluctuations of the Cosmic 
\m) \ Microwave Background (CMB) radiation and the distribution of large scale structures [I]- The basic outcome of 

■ simplest single-field slow-roll inflation models is the statistical isotropy, the Gaussian and (almost) scale-invariant 
O i] power spectrum Q. These statistical properties can be tested by precise measurements of the CMB temperature 
JL . anisotropies. 

^ ' The observational data provided by the Wilkinson Microwave Anisotropy Probe (WMAP) Q showed the evidence 

, of the scale dependence of the power spectrum, whose property has been used to discriminate between a host of 

' inflationary models. Interestingly, the WMAP data also indicated the deviation from the Gaussian perturbations 
and they gave us a hint of the statistical anisotropy According to the recent results of Planck, the Gaussian 
perturbations are still consistent with the data, but the statistical anisotropy remains @. In the light of these new 

^ I results, it is worth investigating a possible statistical anisotropy based on concrete theoretical models. 

■ Naively, the statistical anisotropy implies that vector fields can play an important role during inflation. A mechanism 
' for creating the statistical anomaly at the end of inflation was proposed in Ref. and extended in various ways Q . 
, Moreover, a more concrete model has been proposed in the context of supergravity There, the anisotropic 

r I ' inflation is realized as an attractor and the vector hair survives during inflation. The latter vector hair gives rise 

CO ■ to rich phenomenology [lol-[T2j such as the statistical anisotropy and the cross correlation between the curvature 

' perturbations and the primordial gravitational waves (see Refs. [H, [i| for reviews). In particular, the latter would 

[ imply the correlation between the temperature fluctuation and the B-mode polarization jl5| . 

. Remarkably, subsequent works revealed that vector fields can also induce the large non-Gaussianity [l^. In par- 

^ ' ticular, the non-Gaussianity of curvature perturbations has been investigated in the context of anisotropic inflation 



17|-|22|. which further emphasized the rich phenomenology of the anisotropic inflationary models |23j . As a result 



the Planck constrained the anisotropic inflationary models with vector flelds strongly. In Ref. however, it was 
, pointed out that not only vector flelds but also two- form fields can potentially give rise to anisotropic inflation. In fact. 



it is well known that there are two- form flelds in string theory [2J| . Hence it is natural to explore this possibility from 
the theoretical point of view. One may suspect that there is no statistical anisotropy because the two-form fleld can 
be represented by a pscudo scalar fleld, i.e., axion. However, there remains a possibility that a non-trivial polarization 
of a two-form fleld induces the statistical anisotropy. In this case the anisotropy comes from an expectation value of 
the spatial derivative of the axion fleld. At flrst glance this seems to be odd, but it is a natural framework from the 
picture of two-form flelds. Hence, in this paper, we study this possibility in detail. 

We show that several interesting features are present in our model. Analogous to the results of Ref. @, anisotropic 
inflation can be sustained by the background two-form field. Moreover, as suggested in Ref. [ll| , there exists a prolate 
type anisotropy in the power spectrum of curvature perturbations. In contrast to vector models the non-Gaussianity 
vanishes in the squeezed limit, but the nonlinear estimator /nl of the equilateral and enfolded shapes can be as large 
as the order of 10. Hence, our predictions are consistent with the Planck data and the future analysis may reveal 
these statistical anisotropies. 

The organization of our paper is as follows. In Sec. [Hi we introduce the two-form field model and explain how 
the non-trivial hair remains during inflation. In Sec. Illli we quantize the perturbations of two-form fields and derive 
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their vacuum expectation values on super-Hubble scales. In Sec. IIVI we evaluate various statistical quantities-such 
as n-point correlation functions (n = 2,3,4) of curvature perturbations. Sec. |V]is devoted to conclusions. 



II. ANISOTROPIC BACKGROUND 

We study the background dynamics of anisotropic inflation in the presence of two-form fields. The analysis is similar 
to the case of vector models studied in Ref. Q, but we repeat it for completeness. Here, we emphasize the shape of 
anisotropy is different from that of vector models. 

Let us start with the action 



(1) 



where Mp is the reduced Planck mass, i? is a scalar curvature calculated from a metric g^jy (with a determinant g), 
<j) is an inflaton field with a potential V{(j))^ and the field H^j_^\ has a non-trivial coupling with the inflaton. The 
field i?^i/A is related to a two-form field B^^, as 



(2) 



Without loss of generality, one can take the (y, z) plane in the direction of the two-form field. Then we can express 
B^j^i, in the form 



-B^j^y dx^ A dx^ = v{t) dy Adz , 



(3) 



where v{t) is a function with respect to the cosmic time t. 

Since there exists a rotational symmetry in the (y, z) plane, it is convenient to parameterize the metric as follows: 



ds-" 



Af{tfde + e2"(*) [e-4-(*)da;2 -f- e^^^'\dy^ + dz^) 



(4) 



where a describes the average expansion (the number of e-foldings) and a represents the anisotropy of the Universe. 
Here, the lapse function M is introduced to derive the Hamiltonian constraint. With the above ansatz, the action ^ 
reduces to 
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where an overdot denotes a derivative with respect to t. After taking variations, we set the gauge M — \. 
The equation of motion for the two-form field v is easily solved as 

V = ^/(0)-2e"+4- , 



(5) 



(6) 



where A is a constant of integration. Taking the variation of the action with respect to A/", a, cr, 4> and using the 
solution we obtain the following background equations of motion 



a 



3M2 



a = —3a - 
—3a& 



1 

f 



a 



MI; 



2„-2a+4cr 



-2q+4<t 



(7) 

(8) 

(9) 
(10) 



where the subscript in Vij, and / ^ denotes a derivative with respect to its argument 
Let us introduce the energy density of the two-form field 



Pv 



(11) 
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and the shear T, = a. We also define the expansion rate H = a. First, we need to look at the shear to the expansion 
rate ratio 'S/H, which characterizes the anisotropy of the inflationary Universe. Notice that Eq. ^ reads 



2 Pv 
3M2 



(12) 



In the regime where S becomes negligible, the ratio Y,/H should converge to the value 

— — 

where we used Eq. ([7]) under the slow-roll approximation, i.e.. 



a 



3M2 



(13) 



(14) 



In order for the anisotropy to survive during inflation, must be almost constant. Employing the standard slow- 
roll approximation and assuming that the two- form field is sub-dominant in the inflaton equation of motion (jlOp . 
one can show the coupling function f{(j)) should be proportional to e~" to keep py almost constant. In the slow-roll 



regime, the number of e-foldings a is related to cj), as da 
determined as 



-V{4i)d(f)/{MlV,^). Then the functional form of f{<f)) is 



(15) 



For the polynomial potential V oc 0", for example, we have / = e^""-p . The above case is, in a sense, a critical 
one. What we want to consider is super-critical cases where the energy density of the two-form field increases. For 
simplicity, wc parameterize /(</>) by 



f{<t>) 



(16) 



where c is a constant parameter. Let us consider the super-critical cases c > 1. From the definition (jl6p . we can 
derive the following relation 



V 



f,4> „ 



Then, the condition c > 1 translates into 



f V 



> 1 . 



(17) 



(18) 



Thus, any functional pairs of / and V satisfying the condition (jlSp in some range could produce the two-form hair 
during inflation. 

On using Eq. pT| . the inflation equation pUj) reads 



-Sad) — Vd 



Pv 



ev 1^(0) 



where the slow-roll parameter ey is deflned as 



V 



(19) 



(20) 



In this case, if the two-form held is initially small (pi,/l/(0) ^ ey/c), then the conventional single-fleld slow-roll 
inflation is realized. During this stage / cx e""^" and the energy density of the two-form held grows as p^ oc e2('=~i)". 
Therefore, the two- form field eventually becomes relevant to the infiaton dynamics described by Eq. (|19p . Nevertheless, 
the cosmic acceleration continues because pv/V{(j)) does not exceed ey/c. In fact, if py/V{(j)) exceeds ey/c, the inflaton 
fleld (p does not roh down, which makes pv = A2/((/))~2g-2a+4o-y'2 decrease. Hence the condition py <^ V{(j)) is always 
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satisfied. In this way, there appears an attractor where inflation continues even when the two-form field affects the 
inflaton dynamics [9| . 

Let us make the above statement more precise. Under the slow-roll approximation, the inflaton equation of motion 
(ITUl) reads 



- 3d,^ - ^0 + A2/"^/,0e-2«+4- ~ . (21) 

Using Eqs. ([H]) and (gT]), we obtain 

da a V ^0 

Neglecting the evolution of V, and a, this equation can be integrated to give 

2a-4a+2c f Y rf0 (? V 



where is a constant of integration. Substituting this back into the slow-roll equation (P^ . it follows that 



d(j) _ M^V^ c - 1 M'^V^ r .1-1 



da V c V 



(23) 



(24) 



At the initial stage of inflation {a — — oo), the second term of Eq. (p4)) can be neglected relative to the first term. 
In the asymptotic future {a — )■ oo), the term containing D, disappears. This clearly shows a transition from the 
conventional slow-roll inflationary phase, where 

d^- ~' ^^^> 

to what we refer to as the second inflationary phase, where the two-form fleld is relevant to the inflaton dynamics and 
the inflaton gets 1 /c times slower as 

^ = -1^. (26) 
da c V 

In the second inflationary phase, we can employ Eq. ([23|) with discarding the fl term to rewrite the energy density 
of the two-form field as 

-2c,+4a-2cjj^d4> c-l 

Pv = f = —^evV{(l)) , (27) 

which yields the anisotropy 

S 2c- 1 



H 3 c2 

Moreover, from Eqs. ([7]) and ([8]), the slow-roll parameter e = —H/H^ is related to ey as 



ey . (28) 



a 1 / 1 1 \ 1 

7^0' -oPA= -^v , (29) 



where we neglected the anisotropy and used Eqs. and (PSI) . Thus we arrive at the result 

E 2 c- 1 



H 3 c 



(30) 



Therefore, for a broad class of inflaton potentials and two-form kinetic functions, there exist anisotropic inflationary 
solutions, with "E/H proportional to e. We also confirmed the existence of such anisotropic solutions by integrating 
Eqs. ([7))- ([T(I)l mmrerically. Note that the sign of E/iJ is opposite to that derived for the vector field 
For c = 1, we need a separate treatment. In this case, integration of Eq. ((22|) gives 



2^ T72^(« + "o), (31) 
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where ao is an integration constant. Thus, we obtain the anisotropy 

S _ 1 

Notice that the anisotropy stemmed from the two-form field is a prolate type, while the anisotropy created by the 
vector field was an oblate type. This can be understood by the fact that the vector extending to the x-direction 
speeds down the expansion in that direction, while the two-form field extending in the {y, z) plane speeds down the 
expansion in the (y, z) direction. 

Before entering the study of perturbations, we should note the following important attractor mechanism [25j . Taking 
a look at the result ([^5]) . we find that the relation 

c r — d4> 

f^e-' cx e-" , (33) 

holds during the second anisotropic inflationary phase. Recall that this is the critical behavior. As we will see in the 
next section, this attractor gives rise to the scale-invariant spectrum of the two-form field. 



III. PERTURBATIONS OF TWO-FORM FIELDS 



From the phenomenological point of view the anisotropy of the expansion rate needs to be sufficiently small, so it 
is a good approximation to neglect the effect of the anisotropic expansion. However, we cannot ignore the effect of 
the two-form hair. Actually, in the next section, we will show several interesting results. In this section, we prepare 
some tools for the evaluation of n-point correlation functions of curvature perturbations. 

Let us consider the two-form field given by the action 



I 
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In the above action, there exists a gauge invariancc under the gauge transformation 



(34) 



(35) 



Here, since we have the redundancy + d^X , the parameter ^ can be restricted to be = 0. 

For the derivation of the perturbation equations of the two- form field, we consider the isotropic background described 
by the metric 



where t ~ j a ^dt is the conformal time. From Eq. ^ it is convenient to perform the (3-(-l)-decomposition 

diBjk + djBki + dkBij , 
where a prime represents a derivative with respect to r. The interacting action (1341) reads 



Taking the variation of this action with respect to i?oi, wc obtain 
Variation with respect to By leads to 



d_ 

dr 



Then Eqs. (gOl) and (gT]) reduce to 



9, [B^ + d,Bjo - djB.o] = , 



d_ 



-fi^){Bl^+^,B,o+^,Bo^} 



— /2(0) {d,B,k + djBM + dkB,j} 



= 0. 



(36) 



(37) 
(38) 



(39) 

(40) 

(41) 

(42) 
(43) 
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Since we have the gauge degree of freedom 

SBij = (9,Cj - dj^^ , d'Ci = , 

choosing the parameter 
we can take the gauge 

^^B,, = . 

Similarly, by taking into account the following gauge transformation 



we can set 



From Eq. (j42]) it follows that 



B^o = 0. 



Then, Eq. gS]) reduces to 



d_ 



+ f{^)dlB,,=0. 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 

(50) 



Let us check the degrees of freedom. The components Bij have 3 degrees. There are 2 gauge conditions Bij j ~ 0. 
Hence, we have one degree of freedom. This is the reason why we can map the two-form field to a scalar field. 
However, there is an important difference from the scalar field, namely, there exists a polarization in our case. Using 
the polarization tensor cij = ^^ji with ki^ij = 0, we can expand the anti-symmetric tensor field as 



B^J (r, x) = J [ak x(r, k) e„- e*^'" + c.c] 



(51) 



where ak is the annihilation operator. We take the polarization tensor to be eij — {km/k)€mij, with the normalization 
condition 



e.je'^ =2, 



(52) 



The mode functions x(r, k) is known by solving Eq. (j50p . For convenience we introduce the following variable 

w(T,fc) = ^x(T,fc), (53) 
a 

and parametrize the kinetic function as / = a^. This is always possible during inflation where both (f) and a = e" are 
monotonic functions with respect to time. Then, we obtain 



p{l-p) 



u = 0. 



(54) 



Hence, for p = — 1 or p = 2, we obtain the scale-invariant spectrum for the two-form field. Although either choice 
is possible, we set p = —1 because, as is shown in the previous section, this is an attractor value realized during 
anisotropic inflation. For p = — 1, we can deduce the mode functions as 



Ha 



(1 + ikr) e 



-ikr 



(55) 



and X = a^w. 
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Now, it is convenient to define 

/ / , 

/ / / 

where Ey^ and SEij correspond to the background value and the perturbation of tlie two-form field respectively. Note 
that Hijk can be negligible on super- Hubble scales. We perform the Fourier transformation of the perturbation SEij, 
as 



^^^3= / 77^^'''-^^^'^3{r.k). (58) 



(2^2 

On super-Hubble scales the Fourier modes are given by 

(t, k) = [cik + a^_f}j £k f^ij , (59) 

where 



In the last approximate equality of Eq. ([60]) we used the solution ([55]) in the limit r on the de Sitter background. 
The vacuum expectation value of the field SEij is given by 

The Infrared (IR) modes are characterized by fc^ < k < kf, where ki and kf are the wavenumbers which crossed 
the Hubble radius at the beginning and at the end of inflation respectively. Since the integral ^ = \n{kf/ki) is 
equivalent to the number of e-foldings N — ln{af/ai) on the de-Sitter background, it follows that 



i^Ef,) = ^N. (62) 



On super-Hubble scales the total two-form field is given by 



with the variance (|62p of the perturbation SEij . 



^classical ^ ^ g^^^ ^ (gg) 



IV. STATISTICALLY ANISOTROPIC NON-GAUSSIANITY 



In this section we estimate the statistical properties of our model, in particular, the scalar non-Gaussianity. Through 
this section, the anisotropy is assumed to be sufficiently small. We derive the interacting Hamiltonian by expanding 
the action around the anisotropic background solution. We compute correlation functions according to the in-in 
formalism by neglecting the anisotropic expansion of the Universe. The calculation is analogous to that carried out 
for the vector field in Ref. This prescription can be justified by more rigorous calculations (see e.g., Ref. [13]) • 



A. Power spectrum 



We first calculate the power spectrum of the comoving curvature perturbation C (see Refs. [2a] for its definition). In 
our case the curvature perturbation can be written as the sum of the "unperturbed" field C^"^ and the contribution 
SC coming from the two-form field. We decompose the field C'*^^ into the Fourier components 

C'"' - / (^^"-^Cr^ , = + C^al, , (64) 
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where the annihilation and creation operators satisfy the commutation relation 

At leading order in slow-roll we have the following solution [l^] 

(0) ^ Hjl + ikr) 

The total power spectrum Vq is defined by the two-point correlation function of (, as 

27r2 



(65) 



(66) 



(67) 



where is the Fourier component of (. The power spectrum can be written as the sum of the two contributions ^f"' 
and S(, as 



-Q ~ I- Q ui-Q . (68) 
Using the solution long time after the Hubble radius crossing (t — )• 0), the first term in Eq. (|68p reads 



(0) 

C 



87r2eM2 



(69) 



The next step is to derive the second contribution 5Vc^ — <5(0|CfeiCfc2 10) from the two-form field. The interacting 
Lagrangian following from Eq. is 



12 d(j) 



(70) 



where ( ) represents the background value. Since the function / is given by / = cxp(/ d^/v^Mp), we can deduce 
the following relation 



(71) 



where we used 5(j) = \/2eMpC^''^ ■ Note that there is no distinction between e and cy because we are considering the 
situation c ~ 1. Thus, we obtain 



1 



^int 



(72) 



where, in the second line, we employed the solution / = a^^ and Eqs. ([55]) - ([57]) . 

The interacting Hamiltonian H-mt is related to Lint, as i7i„t = — J (PxLmt = Hi + H2, where Hi and H2 follow 
from the first and second terms of Eq. ([7^ . Substituting Eqs. ([55)1 and into Eq. ([7^ . we obtain 



Hi 

H2 



2E,, 



cPk 5£y,iT,k)C^^liT), 
d^kd^p 



2H*T^ J (27r)3/2 



(73) 
(74) 



Using the in- in formalism [27[ , the two-point correction following from the interacting Hamiltonian Hi can be evaluated 
as 



dri 



dT2 



'^1(^2)1 |o) 



9e2Af4i/4 



n 



dn 



2-K 



6'^^\ki+k2) 
Ki 



{0\6£y,iTi,ki)6£y,iT2,k2)\0) 



47r2e2 Mi 



(75) 



9 



where Oki.x is the angle between ki and the a;-axis. The two integrals in the first line of Eq. ((75)) have been evaluated 
at the super-Hubble regime characterized by — /c^r < 1, from which T^ia,i = — with i = 1,2. In the second line 
of Eq. ([75]) the upper bound ti of the second integral has been replaced by r by dividing the factor 2! because of the 
symmetry of the integrand. Note that we also used the following relation [l^ 

[Cr'M,e^(rO] --^^^g^^^^H^ + fcO, (76) 

which is valid in the super-Hubble regime. One can show that the integral /Ij^y^. d,Ti (t^ — Tf)/Tf is equivalent to 
— ln(— l/fcir) under the approximation — fc^r <^ 1 [l^, where Nk^ is the number of e-foldings before the end of 
inflation at which the modes with the wavenumber ki left the Hubble radius. Since fci = — ^2, we used the notation 

Nk, = iVfc, = Nk. 

Using the power spectrum (p^ and the slow-roll relation 3MpH^ ~ V, the two-form field gives rise to the correction 
to the power spectrum 



(77) 



The interacting Hamiltonian ([71)1 has a contribution to dV^ with Eyz replaced by the IR solution SEij with the 
expectation value (|62p . Taking into account this contribution, the total correction to the power spectrum can be 
derived by replacing Ey^ in Eq. dTT]) for i;<=iassicai defined in Eq. Using the notation E^ = | £;ciassicai | ^ ^^le total 

power spectrum of the curvature perturbation reads 



Vc = ^ (1 + 121 cos' Ok,^,) , where ^ = ^ • (78) 

In contrast to the vector case where the anisotropy is oblate (lol . , we now have the prolate anisotropy. 

The statistics of the WMAP anisotropics uses the parametrization ~ 'P^^\l + cos^ fl^.v), where V is a 
"privileged" direction close to the ecliptic poles [H, ll^- The WMAP data provides the constraint — 0.29 ± 0.031 
[30| . From Eq. ([75)) the anisotropy parameter is 

= 12/iV| , (79) 

from which we obtain 



E 



2 



/ = 2.3 X 10" 



(It) m ■ 



Note that the quantity le = E'^/{2V) characterizes the ratio of the energy densities of the two- form field {E'^/2) and 
inflaton {V), which is much smaller than 1 from Eq. (|5(I)) . 



B. Bispectrum 



The three-point correlation of ^ can be evaluated by using the in- in formalism along the same line of Ref. |i8( | . The 
tree-level contribution coming from the interacting Hamiltonian ()73p is given by 



S{o\CkMks\o) 



dri 



-i/fci 



dT2 



-1/fc. 
3 

lOSeSM^i/e H 

^ 2—1 



-l/k. 

dn / 



^3 3\ 



Ci?Ci°^Ci3'(T),i/2(ri)] ,iJi(r2)] ,H,{t,)\ |0) + 2perm. 

d^p 
(27r)3/2 



X (0|(5£ij(ri, fci - p)5Eij{Ti,p)5Ey:,{T2, k2)S£yz{T3,k3) + 

SE,j{T2,k2 - p)S£^J{T2,p)S8y:,{T3,k3)S£y:,{Tl,kl) + 
SEij{T3,k3 - p)S£^J{T3,p)S£y:,{Tl,kl)S£y:,{T2,k2)\0) 

COS Oki,k2 COS dki,x COS dk^.x 



8V27r3/2e3M6 



5'^^\ki + k2 + k3)Nk,Nk.,Nk. 



2 perm. 



(81) 
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where "2 perm." represents two terms obtained by the permutation. In the last hne of Eq. (|8T|) we used the relation 
eij{ki)tij{k2) = 2cos0fcj^fc2. The loop contribution following from the product of the three interacting Hamiltonians 
H2 provides the bispectrum in which Ey^ of Eq. ([8T|) is replaced by the IR solution 5Eij with the variance ([62|) . Then 
the total anisotropic bispectrum Sc, defined by 5{Q\Ck,Ck2C,kM = B^5'^^\ki + fea + feg), reads 



cos0ki.k^'^oa0ki.xcos0k2,x , „ 

— — h 2 perm. 



(82) 



where / = E^/{2€V) is given by Eq. (|80l). 

We define the non-linear parameter /nl according to the relation 
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by which we have 



/nl = 60/ 2^ 3 3 [ri cos 9k^M2 cos .^^ cos Bk^ ,x + cos 6'fe2 ^fc, cos Ok^ ,x cos a/fc, 

[I C) ^ + ^2 < ''3 

r2Cos6'fc3,fcj cose'fc3.a;Cose'fcj,:c] , (84) 



where 



In the following we employ the approximations {V^^ ~ (7'^"'')^ and Nk^ — — = Nqmb- We also fix r2 = 1 
and define the angle /? = tt — 612 in the range < /? < tt (i.e., < < 2). In this case there is the following relation 

= 2(1 -cos^). (86) 

The local, equilateral, and enfolded shapes correspond to (i) /? — > 0, r3 0, (ii) /3 — tt/3, — 1. and (iii) /? — > tt, 
r3 — >■ 2, respectively. 

Let us consider the situation in which the angle between fci and the x-axis is given by 7. On using Eq. (|80|). the 
non-linear parameter (|84|) reduces to 

/ ^CMB \ 



where 

^^(^3,7) 



2 + 4 



cos (3 cos 7(cos /3 cos 7 + sin f3 sin 7) + - (cos /? cos 7 -f sin /3 sin 7 — cos 7)^ 



From Eq. we can express /3 in terms of r3, as cos/3 = 1 — r|/2 and sin/3 = ?'3-\/l ~ ^1/4, so that /nl is a 
function of for a given value of 7. The non-linear parameters for the local, equilateral, and enfolded shapes are 
given, respectively, by 



= 7.5(1^) (^j/3^sinS, (89) 

.oquil _ o '7 ( 9*\ ( NqMB \ 



/Sr'"" = 29.8(|i-)(^)co.-T, (91) 

where, in the local case, we expanded /nl^' around /3 = 0. The local non-linear parameter ((89|) vanishes in the 
squeezed limit /3 — >■ 0, which is a distinctive feature of our model. The reason why this happens is that, unlike the 
vector models, /nl is proportional to the inner product of two vectors fc.j and kj. In Eq. (1841) the squeezed limit 
corresponds to the case in which the angles 0k2,k3 and Ok^.ki approach tt/2 with 0. 
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FIG. 1: The non-linear estimator /nl versus rs — fcs/fci for a number of different values of cos 7 with g* = 0.1 and A'^cmb = 60. 
The left and right panels show the plots for the angles < 7 < 7r/2 and 7r/2 < 7 < tt, respectively. The local, equilateral, and 
enfolded limits correspond to ra = 0, rs = 1, and rs = 2, respectively. For 7 close to 7r/2, the equilateral non-linear parameter 
is largest. For 7 close to or tt, /nl has a maximum at ra = 2. 



From Eq. ((90|) the equilateral non-linear parameter does not depend on the angle 7. For = 0.3 and A^cmb = 60, 
/nl"' is as large as 10. The enfolded non- linear parameter (|9T|) depends on 7. For cos^ 7 = 1, =0.1 and A^cmb ~ 60, 

jcnfoldcd l^j.gg j^g 

In Fig. [T] we plot the non-linear parameter (j87| versus ra (0 < rs < 2) for g» = 0.1 and A^cmb = 60. The left 
panel and right panel correspond to positive and negative values of cos 7, respectively. In the local limit (ra — 0), 
the estimator /nl vanishes for any value of 7. For the angle 7 close to 7r/2, /nl has a maximum at the equilateral 
configuration (r-i = 1). With the increase of | C0S7I, however, the enfolded estimator gets larger. In particular, for 7 
close to or TT, /nl has a maximum at ra = 2. 



C. Trispectrum 

The four-point correlation function of ^ corresponding to the tree-level contribution is given by 



<5(oiCfciCfc.Cfc3a4io> 



-i/fci 



-l/fc3 



dTA 



-l/ki 



Cl?Cci3^Ci"?(r),i/2(ri) ,if2(r2) ,H,{t,) ,H,{t^) |0) + 5 perm. 



1 E^^ ^ r 



d^pd^p' 
(27r)3 



x(0|(5£j-,(ri, fci ~ p)S£ij{Ti,p)d£ij{T2,k2 - p')5£ij{T2,p')8£yz{Tz,k'i)5£yz{Ti,ki) 
+ 11 perms. 1 0) 

9 



8 -(277)3 g4jv^8 



,3,3,3 cos^fcj^fcj, cos6'fc2,fcj3 cos6';ei,a; cos9k2,x + 11 pcrm. 

"'1"'2"'13 



(92) 



where fc^- = fe^ + fc^. The loop contribution, which follows from the product of the four interacting Hamiltonians iJ2, 
gives rise to the four-point correlation function where Eyz in Eq. (IMI) is replaced by the IR solution 5Eij with the 
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variance ((62)) . Defining the total anisotropic trispectrum Tq, as 5(0|CfciCfc2Cfc3Cfc4 10) = "Tc (^^^-'('"i + ^2 + fcs + ^4), we 
obtain 



Z;>3 Z-«3 i->3 
."'l'>'2"'13 

We introduce the non-Hnear estimator tnl according to the relation 



cos0fc^,fcj3 cos Bk^Ms coB9ki^xCos9k2,x + 11 perm. 



(93) 



Tc = (27r)3(7'c)^^ ( + 11 perm. ) . (94) 



In the squeezed limit characterized by fci2 — >■ 0, the non-linear estimator reduces to 



local 
'NL 



= -lOS/iVl^A^ljcOsflfci^fcj^ COsdk3,ki2 cos 9ki,x cos 6k:,, X + cos 6ki,ki2 C0sdkiM2 COS 9ki,x COS 9ki,x 

+ COS 9k2M2 cos 9k3,ki2 cos 9k2 ,x cos Ok^-x + cos 9k2 M2 cos 9k^,ki2 cos 9k2 ,x cos Bk^ ,x\ , (95) 

where we used the approximation (7-'^)'^ — {J'^^^^Y ■ Since the angles between the vectors ki (i ~ 1,2,3,4) and fci2 
approach 7r/2 for fci2 — )■ 0, the estimator t^^^ vanishes in this limit. 

We also consider the regular tetrahedron limit, i.e., k\ ~ ki^k-^ = k4 = ki2 = fci4 = k (see e.g., figure 2 of Ref. [3l| 
for illustration). For this configuration, the angle between fcia and fci is 7r/4 with ki^ ~ ^/2k. We also focus on the 
case in which the direction of ki is the same as that of the x-axis. Then the trispectrum (|93p reads 

7-^cquii ^ _54(^ ^ l)n^IiV^''^fNl -1 , (96) 
by which the non-linear estimator can be derived as 



r^r' ^ -4.1 X 10^ 



Unlike the local shape, '^^^ ^'^ '^^ ^^^'^ order of 10^-10'^. 



V. CONCLUSIONS 



For the models in which the inflaton field (j) couples to an anti-symmetric tensor Bfj_^, we showed that anisotropic 
inflation occurs for the coupling f{(j)) given by Eq. In this case there is an attractor solution along which the 

ratio of the anisotropic shear E to the Hubble parameter H is proportional to the slow-roll parameter e. Even for the 
super-critical case in which the coupling f{4>) is generalized to Eq. (jl6p with c > 1, there is the regime of anisotropic 
inflation where T,/H is nearly constant with / proportional to a~^. The anisotropy induced by the two-form field 
corresponds to the prolate type (the expansion of the Universe slows down in the {y,z) plane), in contrast to the 
oblate type stemming from the vector field. 

The presence of the two-form field coupled to infiaton gives rise to modifications to statistical quantities observed 
in CMB temperature fluctuations. From the action ([M)) we derived the interacting Hamiltonians ([75)) and ([71)1 
between curvature perturbations and the two-form field. We evaluated the rt-point correlation functions (n = 2,3,4) 
of curvature perturbations by using the in-in formalism of quantum field theory. The 2-point correlation function, 
i.e., the power spectrum, is given by Eq. ([75| . where = 121 parametrizes the strength of anisotropy. Even if the 
energy density of the two-form field is very much smaller than that of infiaton, the parameter can be of the order 
of 0.1, as suggested by the WMAP data 

In Eqs. (j82p and (j84p wc derived the three-point correlation function (bispcctrum) and the non-linear estimator 
/nl, which exhibit a number of interesting properties. By considering the triangle of three momenta (fci -I- fc2 + ^3 = 0) 
with fci = fc2, we showed that /nl can be expressed by functions of ~ fcs/fci and the angle 7 between fci and the 
X-axis. In the local, equilateral, and enfolded limits, the non-linear estimators are simplify given by Eqs. (j89p . (j90p . 
and (PT|) . respectively. We found that f^ff^ vanishes in the squeezed limit (ra — > 0), whereas /^l"' ^iid f^°^^'^'^ can be 
of the order of 10 (see Fig.[T|). These results are consistent with the recent constraints by Planck, i.e., /^l*^' = 2.7±5.8 
and /^r' = -42 ± 75 (68 % CL) 0. 

The four-point correlation function Tq (trispectrum) has been also computed in Eq. (j93p . Defining the non-linear 
estimator tnl as Eq. we found that tnl vanishes in the squeezed limit (fci2 — >■ 0). However, for other shapes 
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such as the regular tetrahedron, |tnl| can be of the order of 10^-10^. This is an interesting property by which our 
scenario can be distinguished from the vector case as well as other models with large non-Gaussianities. 

It will be of interest to understand the physics of a dipole-type anisotropy suggested by the Planck data . Although 
there is a phenomenological description of this type of anisotropy [3^ , no physically well- motivated models are present 
to our best knowledge. Recent attempt to explain the dipole-type anisotropy by a contrived geometrical set up is 
intriguing (33l |. but it still lacks a consistent dynamical picture. Since our framework is natural and consistent, it 
would be great if our mechanism is generalized to explain the origin of the dipole-type anisotropy observed by Planck. 
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